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Abstract. Precise predictions of atomic energy levels require the use of QED, espe- 
cially in highly- charged ions, where the inner electrons have relativistic velocities. We 
present an overview of the two-time Green's function method; this method allows one 
to calculate level shifts in two-electron highly- charged ions by including in principle 
all QED effects, for any set of states (degenerate, quasi-degenerate or isolated). We 
present an evaluation of the contribution of the screened self-energy to a finite-sized 
effective hamiltonian that yields the energy levels through diagonalization. 



1 Experiments and Theory 

Experimental measurements of atomic energy levels provide more and more 
stringent tests of theoretical models; thus, the experimental accuracy of many 
measurements is better than the precision of theoretical calculations: in hydro- 
gen in helium [§J|], and in lithium-like uranium ||] and bismuth |J. The 
current status of many precision tests of Quantum-Electrodynamics in hydrogen 
and helium can be found in this edition. 

Furthermore, highly-charged ions possess electrons that move with a velocity 
which is close to the speed of light. The theoretical study of such systems must 
therefore take into account relativity; moreover, a perturbative treatment of the 
binding to the nucleus (with coupling constant Za) fails in this situation 0. 
Perturbative expansions in Za, however, are useful in different situations (see 
H for a review, and articles in this edition |p|, |To| , pT|Jl^] ) . 



2 Theoretical Methods for Highly-Charged Ions 

There are only a few number of methods that can be used in order to predict 
energy levels for highly- charged ions within the framework of Bound-State Quan- 
tum Electrodynamics j^J: the adiabatic ^-matrix formalism of Gell-Mann, Low 
and Sucher |l4[] , the evolution operator method [ p^5|Jlr^] , the two-time Green's 
function method |l7j and an interesting method recently proposed by Lindgren 
(based on Relativistic Many-Body Perturbation Theory merged with QED) Q. 
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All these methods are based on a study of the some evolution operator or propa- 
gator; the two extreme times of the propagation can be both infinite (Gell-Mann- 
Low-Sucher), one can be finite and the other infinite (Lindgren), and both can 
be finite (Shabaev). 

But among these methods, only two can in principle be used in order to apply 
perturbation theory to quasi- degenerate levels (e.g., the 3 Pi and x Pi levels in 
helium-like ions): the two-time Green's function method and Lindgren's method 
(which is still under development). Both work by constructing a finite-sized 
effective hamiltonian whose eigenvalues give the energy levels |L9j. 

The two-time Green's function method has the advantage of being applicable 
to many atomic physics problems, such as the recombination of an electron with 
an ion [|o), the shape of spectral lines |H| and the effect of nuclear recoil on 
atomic energy levels [p2.23 . 



2.1 Overview of the Two-Time Green's Function Method 

We give in this section a short outline of the two-time Green's function method. 
The basic object of this method |24j represents the probability amplitude for N 
fermions to go from one position to the other, as shown in Fig. H. 




Fig. 1. The 2-particle Green's function is the amplitude for going from one state of 
two particles to another state 



The corresponding mathematical object is a usual iV-particle correlation 
function between two times: 

a[...a> N (*^ • • • ,XiV, t; X l5 . . . , XiV , t ) (1) 

= (r2|T# ai (x 1 ,t)...#^(x iV ,t) 

X# aJv (x , iV ,0--^ai(xi,Ol«> ) ( 2 ) 

where \Q) is the vacuum of the full Bound-State QED Hamiltonian and 
where the quantum field & is defined as the usual canonical electron-positron 
field evolving under the total hamiltonian in the Heisenberg picture ]l3[ . 

A remark can be made here about Lorentz invariance: the above correlation 
function (or propagator) displays only two times, which are associated to many 
different positions. A Lorentz transform of the space-time positions involved 
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therefore yields many different individual times (one for each position); thus, 
the object (|l|) must be defined in a specific reference frame. And this reference 
frame is chosen as nothing more than the Galilean reference frame associated to 
the nucleus, which is physically privileged. 



Fundamental Property of the Green's Function 

The TV-particle Green's function is a function of energy simply defined through 
a Fourier transform of Eq. (Q) : 

a 7V (x 1 ,...,x 7V ;x / 1 ,...,x / 7V ;^GR) 

= 1 JdAt e iEAt 5^(xi , . . . , x N , At; xi , . . . , x' N , t' = 0) . (3) 

This function is interesting because it contains the energy levels predicted by 
Bound-State QED: one can show Q that 

5 7V (x 1 ,...,x iV ;x / 1 ,...,x / iV ;^GR) (4) 

(gfei) • • • ^(x N )\n)(n\t(x f N ) • 
E — (E n — iO) 

(f2\^(x' N ) • - 4(xi)|n)(n|^(xi) • • • Hx N )\f2) 



Eigenstates |n) of H 
with charge — N\e\ 



(-d n2+1 E 



Eigenstates |n) of H 
with charge +iV|e| 



E - (-E n + i0) 



where is the vacuum of the total hamiltonian ^ is the usual second- 
quantized Dirac field in the Schrodinger representation and E n is the energy 
of the eigenstate n of H. The poles in E with a positive real part are exactly 
the energies of the states with charge — iV|e|, which are physically the atomic 
eigenstates of an ion with N orbiting electrons (The charge of the nucleus is not 
counted in the total charge.), as shown graphically in Fig. ||. Such a result is 
similar to the so-called Kallen-Lehmann representation p^| . 

In order to obtain the energy levels contained in (||), we must resort on a per- 
turbative calculation of the correlation function ([l]), which belongs to standard 
textbook knowledge |2(| . The position of the poles of (||) must then be mathe- 
matically found. It is possible to construct an effective, finite-size hamiltonian 
which acts on the atomic state that one is interested in; the eigenvalues of this 
hamiltonian then give the Bound-State QED evaluation of the energy levels JlOfl - 
This hamiltonian is obtained through contour integrations. 



2.2 Second-Order Calculations 



The current state-of-the-art in non-perturbative calculations (in Za) of atomic 
energy levels within Bound-State QED consists in the theoretical evaluation 
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|gN = 2 (E )|^ 



Poles at the energy 
levels of a 2-electron ion 



1s* 



1s, 2s 



Energy 



Fig. 2. The 2-particle Green's function contains information about the atomic energy 
levels of a 2-electron atom or ion 



of the contribution of diagrams with two photons (i.e. of order a 2 , since the 
electron-photon coupling constant is e). For instance, for ions with two electrons, 
the screening of one electron by the other is described by the six diagrams of 
Fig. | 




However, most of the calculations of contributions of order a 2 were, until 
very recently, restricted to the very specific case of the ground-state (see [^7j 
for references). The extension to the calculation of the energy levels of quasi- 
degenerate states represents one of the current trends of the research in the 
domain of non-perturbative (in Za) calculations with QED. 

We have calculated the contribution of the screened self-energy (first and 
fourth diagrams of Fig. ||) to some isolated levels in p7| , p8 .29,30|]. When energy 
levels are quasi-degenerate (e.g., the 3 Pi and x Pi levels in helium-like ions), the 
two-time Green's function method allows one to evaluate the matrix elements of 
the effective hamiltonian between different states; for the first diagram of Fig. ||, 
we obtain the following contribution to this hamiltonian (The two electrons on 
the left are denoted by ni and n2, and the two on the right by n[ and n' 2 , and 
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other notations follow.): 



£(-W-( £ 

P,P' l k*n> pt(1) 



( n P(i) n P(2)\Sl(e np(lv e n ^^)\n f p ^ 1) n f 1 



P'{2)/ 



1 



• £ ( n ^(i) n ^(2)I^U^n P(1) ,e n / p/(i) )|n / p/(1) n / p/(2) )) (5) 

5 ple np (( n P(l)l^b)l n P(l)) 

x(n P(1) n P(2) |/(p-£: n / p/(i) )|n p/(1) n p/(2) )) 

■ 3p/| e ; ((n P(1) n P(2) |/(£ np(1) -P / )I^p/(i)^p/(2)> 
n p'(i) v 

x(n p/(1) |i:(p / )|^p/ ( i))) 



+0^(^-4°))], 

where we made use of standard notations [ p7| : is the energy of the Dirac state 
fc, (— l) pp is the signature of the permutation P o P' of the indices {1,2}, Z 1 
represents the self-energy, and I represents the photon-exchange: 

(ab\I(u)\cd) =e 2 J d 3 x 2 [^ a ( Xl )a^ c ( Xl )] (6) 
x [^l(x 2 )oL u ^ d {x 2 )}D^ v (u\ Xi - x 2 ) 

(a\s{p)\b) eILv <qfc|JMIfcfc) (7) 

Vl Wl 7 2m J ^ e k (l-iO)-(p-u)' KJ 

where a, 6, c and d label Dirac states, and e is the charge of the electron; are 
the Dirac matrices, and i(j denotes a Dirac spinor; the photon propagator D is 
given in the Feynman gauge by: 



exp (i\r\ \/uj 2 — ji 2 + iO^j 



D vv > (u\r)= g vv > — j-j , (8) 

where \i is a small photon mass that eventually tends to zero, and where the 
square root branch is chosen such as to yield a decreasing exponential for large 
real- valued energies u. In Eq. (||), d x \ XQ is the partial derivative with respect to 
x at the point xo, and the skeletons of the screened self-energy diagrams with a 
self-energy on the left and on the right are defined as: 

(np(i)7lp (2 )|^(p,l? , )I^P(l)^P(2)) = 

(np (1) np (2) \I{p-p')\kn' pf{2) ) ^ _ p , (k\S(p')\n' pt(1) ) , 
(np ( i ) np ( 2)|5^(p 5 p / )I^P(i) n P(2)) = 
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(n P(1) \U(p)\k) _ _ ^ (kn P(2) \I(p - p') |n' p(1) n' p(2) ) . 

The terms of order a 2 (E^) — En^) are not included in the above expression 
because they do not contribute to the level shift of order a 2 in which we are 
interested. (They contribute to higher orders, as can be seen in the particular 
case of two levels |3l|, p. 27].) 

This expression is only formal and must be renormalized |27| ; angular in- 
tegrations can then be done and numerical computations can be performed in 
order to yield the Bound-State QED evaluation of the energy shifts. 

For the contribution of the first diagram of Fig. || to have any physical mean- 
ing, it is necessary to calculate it together with the vertex correction (fourth 
diagram of Fig. ||). We have obtained the following contribution to the effective 
hamiltonian for the vertex correction: 

^( _1 ) PP/ Yl (^ n P{2)\I(£n P{1) -£n> p/(1) )\i2n f p , {2) ) 
P,P' iiM 

x J_ r (n P(1) i 2 |/(^)Kin / p/(1) ) 

X 2W U K (1 - iO) ~ (e np(1) ~ v)] [e i2 (1 - iO) - (^ p/(1) - u)] 
+0[o?(E<p-EW)] 

where (711,722) and (^x,n 2 ) still represent the electrons of the two states that 
define the hamiltonian matrix element given here, and where the sum over i\ 
and %2 is over all Dirac states. 

3 Conclusion and Outlook 

We have presented a quick overview of the current status of theoretical predic- 
tions of energy levels in highly-charged ions with Bound-State Quantum Elec- 
trodynamics. We have given a short description of the two-time Green's function 
method, which permits the calculation of an effective hamiltonian that can in 
principle include all QED effects in energy shifts. We have also presented the 
specific contribution of the screened self-energy in the general case (isolated lev- 
els, quasi-degenerate or degenerate levels); the expression obtained can serve as 
a basis for numerical calculations of the corresponding effective hamiltonian. 
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